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The two-dimensional SU(iV) quantum antiferromagnet, a generalization of the quantum Heisen- 
berg model, is investigated by quantum Monte Carlo simulations. The ground state for N < 4 is 
found to be of the Neel type with broken SU(iV) symmetry, whereas it is of the Spin-Peierls type for 
N > 5 with broken lattice translational invariance. No intermediate spin-liquid phase was observed 
in contrast to previous numerical simulations on smaller lattices [Santoro et al., Phys. Rev. Lett. 83 
3065 (1999)]. 



The existence of a short range resonating valence bond 
(RVB) spin liquid Q| is one of the central problems for 
low-dimensional quantum spin systems. An RVB spin 
liquid exhibits a finite gap for spin excitations, has only 
short range order, and does not break any lattice sym- 
metry. The search for RVB spin liquid states was mo- 
tivated by the suggestion that such strongly-correlated 
but quantum disordered states can be turned into a su- 
per conducting state upon doping, which may explain the 
mechanism of the copper-oxide super conductors Q . 

An RVB spin liquid is presumably created by strong 
quantum fluctuations which destroy magnetic ordering. 
The simplest construction of spin liquid states is on lat- 
tices with an even number of spins per unit cell, such as 
two spin ladders ||, bilayerj| or coupled plaquettes |j) . 
There a strong coupling within the unit cells leads to the 
formation of weakly coupled spin singlets. 

This mechanism does not apply to lattices with an 
odd number of spins per unit cell, such as the square 
lattice relevant for the high-T c cuprates. In the square 
lattice Heisenberg antiferromagnet quantum fluctuations 
decrease the Neel order as the spin S is decreased whereas 
the ground state remains ordered even for S = 1/2 [§. 
Stronger quantum fluctuations are thus needed and frus- 
trating interactions have been proposed as one route [Q. 
Since frustrating interactions generally cause a sign prob- 
lem for quantum Monte Carlo simulations, numerical cal- 
culations are usually restricted to small lattices and it is 
hard to draw definitive conclusions. Still, it could be es- 
tablished that the ground state of the antiferromagnetic 
Heisenberg model on a triangular lattice is magnetically 
ordered || A disordered ground state was suggested 
f[of for the fully frustrated model on a kagome lattice, 
but the nature of this state is not clear. The only clear 
evidence for an RVB spin liquid in a frustrated system 
so far is for the hardcore dimcr model on a triangular 
lattice 

The route we follow in this Letter is to increase quan- 
tum fluctuations by considering models with higher sym- 
metry than SU(2) and determine the nature of ground 



state once quantum fluctuations are strong enough to de- 
stroy Neel order. While SU(3) symmetric points occur in 
spin S = 1 Heisenberg antiferromagnets with additional 
biquadratic interactions jyl [u], fil|, SU(4) symmetric 
models have had a long history as special points in cou- 
pled spin-orbital models jl6| . Quantum Monte Carlo sim- 
ulations for the SU(4)-symmetric Kugel-Khomskii model 
]l6| suffer from a negative sign problem in more than 
one dimension |l7| and numerical simulations are thus 
restricted to small lattices Jl8|]. For a related antiferro- 
magnetic SU(4) model in which there is no negative sign 
problem, an RVB spin liquid ground state was suggested 
by previous quantum Monte Carlo simulations Jl9| |, as 
the first example of such a state on a non-frustrated lat- 
tice. These models can be regarded as a generalized 
Heisenberg antiferromagnet, and belong to a family of 
SU(iV) models employing one of the simplest represen- 
tations for the symmetry group of the model. In this 
sense they are "minimal" models for physical systems 
with higher symmetries than SU(2), and the results are 
expected to be of relevance as reference models in a wide 
array of applications. 

An SU(A r ) model is specified by an irreducible repre- 
sentation of the operators, i.e., the Young tableau for the 
SU(iV) algebra. Here we consider the series of models 
with a single-row Young tableau. In one dimension the 
ground state of these models is dimerized for any fixed 
number of columns n c of the Young tableau as N — * oo 
pof . In two dimensions, in a bosonic representation 
equivalent to a single-row Young tableau, it was found 
that for large n c and N, the ground state is Neel state if 
N < nn c , where k is a numerical constant |2l] . Read and 
Sachdev p2] | confirmed this result and, by examining the 
continuous version of the model in the large- N limit, ar- 
gued that the ground state becomes a spin-Pcicrls state 
as soon as N exceeds the boundary N* = nn c . Although 
their approach is less reliable for small N, they conjec- 
tured a direct transition with no spin liquid intermediate 
phase even for small values of N (i.e., for small n c ) - 
inconsistent with the numerical observation of a spin liq- 



uid phase in the SU(4) model [|19|. In this Letter, using a 
new cluster quantum Monte Carlo algorithm, we present 
conclusive numerical evidence for a direct transition from 
the Neel state to the spin-Peierls state. No RVB spin liq- 
uid ground state is realized for any N, in contrast to the 
previous simulations on small lattices 

Our SU(iV)-invariant generalization of the Heisenberg 
model pOL E3] can be formally written as 



H= J2 Hrr- = ^ £ ^(r)jf(r'), 
(rr 1 ) <rr') 



(1) 



where {rr') runs over all nearest-neighbor pairs on a 
square lattice, and repeated indices a, (3 — 1, ...,N 
are to be summed over. Throughout this letter, we 
set J = 1 as the unit of the energy. The symbols 
Jg (r) denote the generators of SU(iV) algebra, that sat- 

isfy lJ$(r),JZ(r')] = 5 r ,r> (s^{r) - Sp°(rj) . To 
uniquely specify the model, we have to choose the rep- 
resentation of the algebra. The model examined in this 
letter is the "antiferromagnetic" SU(iV) model, in which 
an operator Jo is represented by a N dimensional matrix 
with the fundamental representation (i.e., a single-box 
Young tableau) on one sublattice and with the conjugate 
representation (i.e., the single-column Young tableau 
with N — 1 boxes) on the other sublattice. We note that 
the Kugel-Khomskii model jl6| uses the same representa- 
tion on both sublattices and that for N = 2 both models 
reduce to the ordinary SU(2) antiferromagnetic Heisen- 
berg model. 

In this representation, the model can be conveniently 
expressed in terms of SU(2) spins with 5 = (N — l)/2. 
The matrix elements of the pair Hamiltonian Hrr> are 
explicitly given by 



(a',p'\Hrr'\a,!3) 



J 

N 



(2) 



where \a, (3) (a, (3 — — 5, —S + 1, • • • , 5) is the simul- 
taneous eigenstates of the z components of SU(2) spin 
operators, S z {r) and S z (r'). We probe for two types of 
long-range order. The first is a generalized Neel state 
with broken SU(iV) symmetry characterized by a finite 
staggered magnetization 



£(-l) r S*(r) 



(3) 



The second is a dimerized state with broken translational 
invariance but no broken spin rotation symmetry, char- 
acterized by an order parameter such as 



D, 



fc =J>-^S z (r)^(r + e;c ) 



(4) 



where e x is the lattice unit vector in the x direction. It 
was argued Q that k = (tt, 0) or (0, tt) is preferred when 
the lattice translational symmetry is broken. 
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FIG. 1: One of the three-fold horizontal graphs for N = 4 
corresponding to A^ rizontal (r,r'), with one of its matching 
spin-configurations. (A different choice of the graph does not 
make any difference because of the symmetrization in Eq. (||) .) 
Open circles denote up-spins and filled circles down-spins. A 
spin-configuration matches a graph if and only if any two 
connected spins are antiparallel. 



Recent developments in quantum Monte Carlo algo- 
rithms [ p3| allow us to perform quantum Monte Carlo 
(QMC) simulations on larger systems and for a wider 
range of N than was possible in previous calculations 
p9| . By splitting each original spin operator into 25 
Pauli spins with 5 — 1/2 |24}| . the Hamiltonian of the 
SU(iV) model in the new extended Hilbert space is ex- 
pressed jl3j as 



J 

'N 



^Horizontal V' ) ' 



(5) 



where the symbol [•••]« denotes the symmetrization with 
respect to the 25 Pauli spins. The symbol A^T 1 ^ ntal 
denotes the operator whose matrix element is one if the 
(N — l)-fold horizontal graph matches the initial and the 
final spin states, while it is otherwise (See Fig. 0). This 
formulation yields the following quantum Monte Carlo 
algorithm: (i) for a given world-line configuration, dis- 
tribute the (N — l)-fold horizontal graphs with the den- 
sity J/N, (ii) construct loops by following world lines and 
horizontal graphs, and (iii) flip each loop independently 
with probability 1/2. The majority of results reported 
in this Letter were obtained with this algorithm. Bet- 
ter algorithms could be constructed recently using the 
framework of the stochastic series expansion [^5[ and the 
idea of coarse-graining p6[ , or directly working with the 
weight equation |]27| in the conventional framework of 
loop algorithms. This latter algorithm was used in a 
part of the calculations. The details of these algorithms 
will be reported elsewhere [ p8| . 

Simulations have been performed at low enough tem- 
peratures to be effectively in the ground state; careful 
checks have been performed for each system size and 
value of N by varying the temperature. We explored 
system sizes L up to L = 128 for N — 3 and 4, and 
L = 64 for N = 5, 6, 7 and 8. The number of Monte 
Carlo sweeps in a typical set of simulation is 3.4 x 10 5 for 
the most time consuming case [L = 128 and N = 4). 

In Fig. H we show the spin structure factor, 5(> ^ = 
£~ d ((Af s ) 2 ), divided by L 2 , which in the SU(iV) language 
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FIG. 2: 



Static structure factors Sr n ,Tt) f° r 2 < TV < 8. The 
straight line representing the power law, S^^L~ 2 oc L~ 2 , 
is drawn for comparison. Estimated statistical errors are not 
shown because they are equal to or smaller than the symbol 
size. The inset presents the data for TV = 3, 4, 5, 6 in the 
linear scale, together with the best fitting curves obtained by 
the method of least-squares. 
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phase, S( n ^)/L 2 converges to the square of the staggered 
magnetization per spin, m s , as L — > oc, while it decreases 
asymptotically to zero, being proportional to L~ 2 , in the 
absence of Neel order. Our results in Fig. || show clear 
evidence for Neel order for TV < 4 and absence of Neel 
order for TV > 5, indicating that the phase boundary of 
the Neel ground state lies between TV = 4 and 5. 

For TV < 4 the staggered magnetization in the ther- 
modynamic limit is calculated from a finite size extrap- 
olation using a second order polynomial in 1/L. A least- 
squares fit based on the data for L > 8, gives 



In the Neel 



m s = 0.1682(6) 
m s = 0.091(3) 



for TV = 3, 
for TV = 4 



(6) 
(7) 



where the numbers in the parentheses indicate the esti- 
mated statistical error (of one standard deviation) in the 
last digit. For TV = 5 and TV = 6, the estimated m 2 agrees 
with zero within the statistical error (m 2 s = —0.0003(6) 
for TV = 5 and m 2 = 0.0004(5) for TV = 6). 

The existence of Neel order at TV = 4 is further con- 
firmed by the correlation ratio, which eliminates contri- 
butions from short-range correlations |29[| . While the 
structure factor St^^) is a sum of two-point correlation 
functions with all distances and therefore contains short- 
range correlations, the quantity 

r Ms (r) = (M a (r)M a (0)> = {e rH M s e- TH M s ) 

does not if r is sufficiently large. We measure this quan- 
tity at a fixed aspect ratio, 0/L, and compute the ratio 
(a'P \L,f3) /Tm b (a/3 | L, fi) with a < a' for various sys- 
tem sizes. As L is increased, this ratio converges to 1 if 
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FIG. 3: The ratio of dynamic correlations of staggered mag- 
netization at two imaginary-time intervals, (3/2 and /?/4, for 
the SU(4) model. As L increases, the value approaches 1, 
regardless of the aspect ratio (3/L. 



the system has long range order. If, on the other hand, 
the system is disordered with a finite gap and therefore 
has a finite correlation length in the imaginary-time di- 
rection, it will converge to 0. Hence this ratio serves as 
a good indicator of the existence of a long range order, 
similar to the well-known Binder parameter. In Fig. ||, 
we plot the correlation ratio for various aspect ratios 
(3/L = 1.0, 1.5 and 2.0 for TV = 4. We can see that 
the correlation ratio approaches 1 and thus establish long 
range Neel order. 

The disagreement of our conclusion (Neel order for 
TV = 4) with that of Ref. |T^] is not in the raw numerical 
data. Their estimates of S^^) agree with ours shown in 
Fig. U within the statistical errors. The disagreement is 
solely due to the small system sizes studied in Ref. [[19). 
The convergence of the magnetization to a finite value in 
Fig. can be seen only for system sizes L > 32, whereas 
the previous simulations were limited to L — 12 

In order to fully answer the question whether an inter- 
mediate spin-liquid phase exists, we next determine at 
which value of TV the ground state starts being dimer- 
ized. Figure f| shows the dimer structure factor, Sj^ = 
L- d {D_ k D k ), for k = (tt, 0) divided by L 2 at TV = 4, 5 
and 6. In the thermodynamic limit, S k /L d should con- 
verge to the squared dimerization per spin. In Fig. |J we 
see a clear power law decay of Sj^/L 2 following L~ 2 and 
thus absence of dimerization for TV ~ 4, but a slower 
decay and up- ward bending for TV = 5 and 6. This sug- 
gests convergence to finite values, although the examined 
systems are not large enough to cover the region where 
iSP Q ^ shows no size dependence. For a more systematic 
analysis, we once more perform a least-squares fit of the 
data of iSjP j for L > 8, using a second order polynomial 
in 1/L. Our results for the spontaneous dimerization per 
site are 0.103(3) for TV = 5 and 0.18(5) for TV = 6 (See 
the inset of Fig. |j). For TV = 4, the same analysis yields 
(£> (7r 0) /X d ) 2 = -0.00002(2), consistent with absence of 
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FIG. 4: The k = (ir, 0) dimer structure factors j for 
N = 4, 5, 6 in logarithmic scale. The inset is the linear-scale 
plot. The solid lines in the inset are the best fitting curves of 
least-squares based on the L > 8 data. 



dimerization 
We also 



compute SP 



at k = (tt,tt) and find 



SP^ ^/L 2 oc L~ 2 , for N = 4, 5 and 6, indicating absence 
of dimerization at this wave vector, consistent with the 
previous suggestion (22). 

In conclusion, our high-accuracy QMC simulations us- 
ing new loop-type QMC algorithms have shown that the 
ground state of the SU(iV) square lattice antiferromagnet 
is the Neel state for N < 4, whereas it is the dimerized or 
Spin-Peierls state for N > 5. No intermediate spin- liquid 
has been observed, consistent with analytical arguments 
p2| but inconsistent with previous numerical simulations 
on smaller lattices pjfl . 

It is interesting to compare the present result with the 
analytical estimate of the phase boundary n c ~ 

0.197V* where n c is the number of the columns of the 
Young tableau. This estimate is supposed to be accurate 
for large N, and is not necessarily justified for the present 
case where n c = 1 but is still surprisingly accurate. For 
our model it would indicate Neel order up to N = 5, 
while we find Neel order only up to N = 4. 

Concerning the existence of true RVB spin liquid states 
without any broken symmetry, our results are essentially 
negative. By showing that the proposed spin liquid state 
in an SU(4) model exhibits Neel order, we are left with 
only the hardcore dimer model on a triangular lattice ]Tl| ] 
as a model with clearly established spin liquid ground 
state |n], ^2|. Since the hardcore dimer model on the 
square lattice does not show any gapped spin liquid phase 
p0| , and in numerical simulations only symmetry broken 
phases were found so far for frustrated lattices, we are led 
to conjecture that a spin liquid state without any broken 
symmetry seems to be impossible to obtain on a bipartite 
lattice with an odd number of spins in the unit cell. 
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